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ABSTRACT 

We use the new cosmological recombination code, CosmoRec, for parameter estimation 
in the context of (future) precise measurements of the CMB temperature and polarization 
anisotropics. We address the question of how previously neglected physical processes in the 
recombination model of Recfast affect the determination of key cosmological parameters, for 
the first time performing a model-by-model computation of the recombination problem. In 
particular we ask how the biases depend on different combinations of parameters, e.g. when 
varying the helium abundance or the effective number of neutrino species in addition to the 
standard six parameters. We also forecast how important the recombination corrections are for 
a combined Planck , ACTPol and SPTpol data analysis. Furthermore, we ask which recombi- 
nation corrections are really crucial for CMB parameter estimation, and whether an approach 
based on a redshift-dependent correction function to Recfast is sufficient in this context. 

Key words: Cosmic Microwave Background: cosmological recombination, temperature 
anisotropics, radiative transfer 



1 INTRODUCTION 

The Planck SurveyoiQ is currently observing the temperature 
and polarization anisotropies of the cosmic microwave background 
(CMB) with unprecedented accuracy, constantly producing new ex- 
citing results. By now it has successfully completed two full sky 
scans, and is about to finished its third. On January 1 1th, 201 1, the 
first release of Planck data became available to the public, e.g. the 
Early Release Compact Source Catalog ( Planck Collaboratio n et al.| 
|201 la[ >, which also contains the Early Sunyaev-Zel'dovich Cluster 
Sample ( jPlanck C ollaborat ion et al.|201 lb} . 

Using Planck data, cosmologists will be able to determine the 
key cosmological parameters with extremely high precision, mak- 
ing it possible to distinguish between various models of inflation 
(e.g. see Komats u et al.|20T0| > by measuring the precise value of 
the spectral index of scalar perturbations, n s , and constraining its 
possible running, n run . In the near future SPTpoip|(McMaho n et al. 
2009 )and ACTPoip1( Niemack et al.|2010| > will provide additional 
small scale ii-mode polarization data, complementing the polariza- 
tion power spectra obtained with Planck . 

These encouraging observational prospects have motivated 
various independent groups (e.g. see Dubrovich & Grachev|2005| 
Chluba & Sunyaev|[2006) |Kholupenko & Ivanchik|2006[ |Switzer| 
& Hirata||2008[ |Wong & Scott||2007[ |Rubino-Martm et al.||2008| 



Karshenboim & Ivanov||2008[ |Hirata||2008t |Chluba & Sunyaevl 
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2008 ; Jentschura 2009 , Labzowsky et al. 2009 ; Grin & Hirata 2010; 
Ali-Haimoud & Hirata 2010b) to assess how uncertainties in the 
theoretical treatment of the cosmological recombination process 
could affect the science return of Planck and future CMB exper- 
iments. It was shown ( Rubino-Martm et al. 2010) that for the stan- 
dard six parameter cosmology in particular our ability to measure 
the precise value of n$ and the baryon content of our Universe will 
be compromised if modifications to the recombination model of 
Recfast ( |Seager et al.]T999l[2000^ are neglected. 

These efforts are now coming to an end, and currently it ap- 
pears that all important corrections to the standard recombination 
scenario have been identified (e.g. see Fendt et al. 2009 ; R ubifio-| 
|Martfn et al.|2010] for overview). This has lead to the development 
of two new independent recombination codes, CosmoRec ( Chluba 
|& Thomas|2"0T0) and HYREc lATTHaimoud & Hirata| l |2010a) , both 
of which include all important corrections to the recombination 
problem, superseding the physical model of Recfast. CosmoRec 
and HyRec allow fast and accurate computations of the ionization 
history on a model-by-model basis, without the necessity of fudg- 
ing. A detailed code-comparison is currently in preparation, how- 
ever, 'out-of-the-box' comparisons indicate agreement at a level of 
~ 0.1% - 0.2% during hydrogen recombination (z ~ 1 100). 

The differences in the free electron fraction with respect to the 
original version of Recfast ( Sea ger et al.|1999| reach the level of a 
ANJN e ~ 1% - 3% (see Fig.[T]for more details). These modifica- 
tions affect the shape of the Thomson visibility function jSunyaev] 
|& Zeld ovich 1970), slightly shifting its position and changing its 
width. This leads to important differences in the theoretical pre- 
dictions of the CMB temperature and polarisation power spectra, 
reaching AQ/C; ~ -4% at / ~ 3000 with respect to the result 
obtained with the ionisation history given by Recfast (see Fig. [2J. 
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Neglecting these corrections leads to biases in the inferred param- 
eters, which are important at the level of a few standard deviations 
for Planck (Rubino-Martm et al.|2010| >. 

In this paper we address the important question about how the 
aforementioned modifications to the ionization history propagate 
into different cosmological parameters. We use the recombination 
code CosmoRec in connection with CAMB (Lewis et al. 2000| and 
CosmoMC ( |Lewis & Bridle 20021 for our computations. 

With CosmoRec it has at last become possible to perform 
model-by-model parameter estimations, including all important 
physical corrections to the recombination process in a detailed re- 
combination calculation. The study by Rubino-M artm et aL|(2010) 
was partially reliant on the results of the multi-dimensional inter- 
polation scheme Rico (Fendt et al. 2009), and the final parameter 
estimations were performed using the newly introduced correction 
function (Rubino-Martm et al. 2010). Also, at that time the recom- 
bination code used for hydrogen ( [Chluba et al.|2007[ l was limited to 
100 shells only, while it is now possible to use effective rates (Ali 



Hami oud & Hirata|2010b| > for up to 500 shells (Chluba & T homas 



2010 1, capturing the correct behaviour in the freeze-out tail of re- 



combination (z < 800). 

Furthermore, we show in more detail how different combina- 
tions of parameters affect the associated biases. In particular, in ad- 
dition to the standard six parameters we allow for a variation of 
the helium abundance, Y p , and the effective number of neutrino 
species, N v , finding that the main biases are strongly reshuffled 
(see Table [2] and [5J. We also extend the analysis to forecast pos- 
sible biases to combined Planck plus ACTPol and SPTpol data 
sets. Finally, we address the question of which physical processes 
are really important in the recombination correction, and whether 
an approach based on a correction function ( Rubino-Mart m et aE| 
2010) is sufficient for future parameter estimation. 

The paper is structured as follows: in Section [2] we describe 
CosmoRec and some of the important corrections to the recombi- 
nation process; in Section[3]we explore the biases produced by dif- 
ferent recombination calculations and their implications for cosmo- 
logical data analysis; finally in Section [4] we discuss the necessary 
accuracy of calculation required for unbiased analysis of the CMB, 
providing some details on a Recfast correction function approach. 



2 RECOMBINATION PHYSICS 

In this section we briefly discuss the different recombination cor- 
rections included in the computations carried out with CosmoRec^] 
For a more general overview on recombination physics we refer the 

(20091 ), |Sunyaev & Chluba) p009) 



interested reader to [Fendt 
and Rubino-Martm et al 



Itetal.l f; 
H2010) . 



We study the precision of CosmoRec using different combina- 
tions of the recombination physics, as well as various settings for 
the solver parameters. The main purpose is to demonstrate that us- 
ing CosmoRec with default setting is sufficient for precise analysis 
of Planck data, and even for a combination of Planck plus ACT- 
pol and SPTpol. We also address this question again in Section [3] 
using CosmoMC. For all results presented in this Section we used 
the fiducial cosmological model given in Table|2] 



CosmoRec is available at http : //www . Chluba . de/CosmoRec 



2.1 Corrections to the cosmological ionization history and 
the CMB power spectra 

CosmoRec allows the incorporation of all known, important, cor- 
rections to the cosmological recombination problem. However, its 
runtime varies rather strongly with the level of detail in the recom- 
bination model, and for parameter estimations using CosmoMC it is 
important to reduce the runtime as far as possible, when using the 
full recombination calculation. It is therefore useful to check which 
corrections need to be accounted for to obtain unbiased results in 
the final parameter estimation. 

In TableQ]we define a set of recombination models to illustrate 
the possible differences in the ionization history. The corresponding 
corrections to the ionization history as a function of redshift are 
presented in Fig.[T] and Fig. [2] shows the associated modifications 
in the CMB power spectra. For the curves given in the left panel 
of Fig. [T] the output from our C++ version of Recfas - ^] was used 
as reference case, while in the right panel the relative differences 
between the alternative models are illustrated, with the CosmoRec 
'default' case defining the reference. 

The left panel of Fig. [T] shows, that the recombination correc- 
tions with respect to the original version of Recfast reach the level 
of a few percent, both during hydrogen and helium recombination. 
During helium recombination (z ~ 1700 - 2200) the largest correc- 
tion ( 2%) is because of the increase in the photon escape rate 

of the 2'P - l'S singlet line mediated by the absorption of helium 
photons in the H 1 continuum ( Kholupenko et al. 2007] |Switzer &| 
|Hirata 2008 , Rubino-Martm et a l.|2008| l. During hydrogen recom- 
bination (z ~ 800 - 1500) details in the radiative transfer of the 
Lyman-series are crucial, including two-photon corrections, reso- 
nance scattering, and Raman-events (see |Chluba & Thomas|2010[ 
|Ali-Hairnoud & Hirata 2010b and references therein for more de- 
tails). The modification in the CMB power spectra related to the 
CosmoRec 'default' model are shown in the left panel of Fig. [2] 

As the right panel of Fig. [T] indicates, the changes in the ion- 
ization history between the different models of Table[T]are already 
very small. For all cases shown the differences in the hydrogen re- 
combination history are AN e /N c < 0.1%. The largest difference 
appears when switching off the diffusion correction to the escape 
probability of the Hei 2'P - l'S resonance 1 Rubino-Martm et al. 
2008), resulting in ANJN t ~ 0.3% uncertainty at z ~ 1800, a cor- 
rection that currently is not included by HyRec. During hydrogen 
recombination, higher level two-photon decays (n > 3) do still lead 
to some ~ 0.1% uncertainty, which appears to be dominated by 
the 4s- Is and 4d-ls process. However, when computing the corre- 
sponding changes in the CMB temperature and polarization power 
spectra (right panel of Fig [2}, it becomes clear that the small dif- 
ference with respect to the CosmoRec 'default' case will not have 
a major impact on the cosmological parameter constraints for cos- 
mic variance limited experiments at / < 3000. The agreement in 
the prediction for the CMB power spectra is better than 0. 1 % for 
all considered recombination models. This error is below the 3/1 
benchmark suggested by Seljak et al. ( 2003]), and as we will see 
below, indeed there is no significant bias introduced when choos- 
ing between the different recombination models defined in Table[T] 

These findings suggests that one can use CosmoRec with the 



5 Recfast++ is part of CosmoRec and can also be downloaded separately 
at http://www.Chluba.de/CosmoRec It reproduces the Recfast result 
I Seager et al. 1999 1, but avoids any switches in the ODE system (see |Fendt[ 
|et al.|2009| for details). Running Recfast++ is similar to running Recfast 
vl.5 jwong et al.|2008| with all hydrogen and helium flags set to zero. 
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CosmoRcc vs Rccfast++ (Rccfast++ is reference) 



CosmoRec different sellings (default case is reference) 




Figure 1. Modifications to the cosmological ionization history. - Left panel: comparison of different recombination models (see Table [TJ computed using 
CosmoRec with Recfast++. - Right panel: relative difference between the different recombination models. Here the CosmoRec default case is used as reference. 



CosmoRec (default) vs Recfast++ 
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CosmoRec different settings (default ease is reference) 
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Figure 2. Changes in the CMB temperature and polarization power spectra. - Left panel: correction for the CosmoRec default case with respect to the standard 
Recfast++ case. - Right panel: additional corrections for different recombination models. Here the power spectra obtained using the ionization history of the 
CosmoRec default case define the reference. In both panels we show the 3/1 benchmark I Seljak et al. 2003) for comparison (dotted line). 



'default' setting to perform accurate cosmological parameter es- 
timations with CosmoMC. Most importantly, for this setting Cos- 
moRec runs in 1.3 seconds per cosmology, and was already tested 
in a wide range of cosmological parameters to assure stability of 
the recombination code. In contrast to the most demanding setup 
('extreme' case of Table[TJ CosmoRec runs about 340 times faster 
in this mode. We also confirmed the precision of CosmoRec by 
comparing directly with the most detailed computation carried out 
using a more elaborate multi-level hydrogen-helium recombination 
code jChluba & T homas 2010). We found differences no larger than 
ANJN C ~ 0.01% at all redshifts. 

We will address the question about the precision of the recom- 
bination model more formally in Section [4] however, we find that 
for accurate parameter estimation using future data from Planck 
, ACTPol and SPTpol the CosmoRec 'default' setting is indeed 
sufficient. However, we would like to point out that a final cross- 



validation of the CosmoRec outputs with independent recombina- 
tion codes ( |Switzer & Hirata||2008| |Hirata||2008| |Grin & Hirata| 
|20 1 0| [Ali-H aimoud et al.|2010[ > will be very important. Neverthe- 
less, we do not expect our conclusions to change very much. 



2.2 Modifications to the freeze-out tail of recombination 



Earlier computations of hydrogen recombination (Chluba et al. 
|2007| > were limited to atomic models with ~ 100 shells. However, 
it was later shown that, as expected ( Chluba et al. 2007 1, the freeze- 
out tail of recombination (z ~ 800) is still affected significantly un- 
til about 300-400 shells are included ( |Grin & Hkata|2010||C^luba 
et al.|2010| >. With the effective rate method ( Ali-Haimoud & Hirata 



2010b I it has now become possible to account for this correction in 
the computation of the CMB power spectra in a very fast way. 

When going from 100 shells to 400 shells, the correction in the 
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Table 1. Different settings for CosmoRec. In all cases the corrections to the Hi 2s-ls two-photon channel were switched on I Chluba & Sunyaev 2006 
Kholupenko & Ivanchik|2006|. Furthermore, we used the effective rates for our hydrogen and helium models with = 500 and n™ 1 = 30 (Chluba 



1 



— — — „, — „ ..„„ ...... ,. eff „„„ ,. cff l— 

Sunyaev|2010||Chluba et al.|2010j. These were computed with the method of Ali-Haimoud & Hirata ( 2010b I. The Hei 2 3 P — l'S intercombination line 
was always switched on (Dubrovich & Grachev|2005} and the effect of H i continuum absorption on the 2'P — 1 1 S singlet and 2 3 P - 1 1 S triplet resonance 
I Kholupenko et al. 2007 , Switzer & Hirata 2008 1 was included using the no redistribution approximation I Rubino-Martm et al. 2008 1. In all cases, the detailed 
history was solved starting at z = 3000 and ending at z = 50 with 500 intermediate points. The solution was completed until z = using the simple Recfast 
ODE system, unless stated differently. Depending on the corresponding settings, the feedback of helium photons was treated according to Chluba & Sunyaev 
(20 10) , including n'D - 1 1 S quadrupole lines and n 3 P — 1 1 S intercombination lines from levels n > 2. The diffusion correction to the escape probability of 
theHEi2'P- l'S resonance was included using the tabulated correction function given in Rubino-Martm et al. | j2008) . Higher level two-photon decays and 
Raman-scatterings were incorporated according to Chluba & Thomas ( 2010), giving results that are in excellent agreement with Hirata 1 2008 1. 



parameter 


default 


HyREC-case^ 


HI transfer 


full 


extreme 


resolved states 


2s-3s, 2p-3p, 3d 


2s-4s, 2p-4p, 3d-4d 


2s-8s, 2p-10p, 3d-8d 


2s-8s, 2p-10p, 3d-8d 


2s-8s, 2p-10p, 3d-8d 


# resolved states 


5 


8 


22 


22 


22 




3 


4 


8 


8 


8 


^Raman 


2 


3 


7 


7 


7 


„HeI 
"max 


2 


2 


2 


5 


5 


Hei diffusion correction 


on 


off 


on 


on 


on 


Hei feedback 


off 


Hf = 2 


off 


Hf = 5 


«f = 5 


Az (PDE-solver) 


20 


20 


2 


2 


0.2 


average runtime 


1.3 sec 


2.2 sec 


38 sec 


50 sec 


350 sec 



This case is intended to reproduce the result of HyRec jAli-Haimoud & Hirata|2010a) using CosmoRec. 



freeze-out tail of hydrogen changes from about AN C /N C ~ 2.8% 
to ~ 1.6% (see for example Fig. 2 in |Chluba & Thomas||2010] >. 
This is expected to lead to some small change in the CMB power 
spectra, since the total optical depth to the last scattering surface is 
slightly modified. However, the difference is very small, and as we 
show here in detail (Section [4]l, for precise parameter estimation 
with Planck it would be sufficient to include only 100 shells to 
the computation. Although this was already suggested by |Rubifk>| 
Martin et al. (2010), here we explicitly show this using CosmoRec, 
however, as the effective rate coefficients can be easily computed 
even up to 500 shells, this does not lead to any additional obstacle 
in the recombination calculation. 

The main implication of this is that changes at the level of 
~ 1 % in the freeze-out tail of recombination are not constrainable 
with Planck. This further suggests that modifications caused by 
collisional processes < |Chluba et al.|2007|[20T0) at z < 800 should 
not matter very much. It was already shown that collisions lead to 
a small acceleration of recombination at z < 800. Since collisional 
rates are very uncertain, this could imply modifications > 0.1% 
( |Chluba et aT1[20Tot . However, since changes ~ 1% at low red- 
shifts (z < 800) do not seem to affect the CMB power spectra at a 
significant level (see Section |4](, the above statement appears rea- 
sonable. Nevertheless, it will be important to check this with refined 
computations of collisional rate coefficients. 



CosmoRec (default) vs Recfast++ (Recfast++ is reference) 
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Figure 3. Modifications to the cosmological ionization history at low red- 
shifts. We compare the output of CosmoRec for default setting with Rec- 
fast++. The difference in the models is only due to the low-z extrapola- 
tion, which is done using the simple Recfast++ ODE system at z < Zextra- 
The models shown here lead to changes in the CMB power spectra with 
AC//C/ < 10" 5 at / < 3000. 



2.3 Computation at very low (z < 200) redshifts 

In the current version of CosmoRec the recombination problem is 
solved until some lower redshift, z cxtra . Below this redshift the sim- 
ple ODE system of Recfast++ is used to complete the solution until 
z = (Chluba & Thomas 2010). The derivatives of CosmoRec at 
Zextra are used to re-scale the derivatives of Recfast++ accordingly. 
Clearly this procedure is expected to introduce some small error to 
the ionization history, however, as Fig. [3] demonstrates the differ- 
ences are small, when varying the value of z cxtra . In particular, we 



found that the associated uncertainty in the CMB power spectra for 
all shown cases is AC//C; < 10~ 5 at / < 3000. For CosmoRec we 
shall therefore use z exlra = 50 for all computations, without intro- 
ducing any significant change to the CMB power spectra. 

This result also suggests that any modifications introduced by 
details in the primordial chemistry {Standi et al.|1996|[T998l|Schle-| 
|icher et al.|2 008) should not have a major impact on the predictions 
for the CMB power spectra, as they are expected to be similarly 
small. Furthermore, at z i 10-20 the process of cosmological 
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reionisation (Barkana & Loeb 2001) after the formation of first 
stars in the Universe (e.g. see |Tegmark et al.p 997; Abe l et al"| 
[T997l[2000l|Yoshida et al.|2007| > is expected to occur. This leads 
to a much larger ambiguity (e.g. see Mortonson & Hu 2008] for 
model-independent estimation) in the cosmological ionization his- 
tory than inherent to detailed recombination calculations. In addi- 
tion, possible changes caused by energy injection from decaying or 
annihilating relic particles (e.g. se ejChen & Kam ionkowski j2004| 
Padmanabhan & Finkbeiner|2005||Zhang et al.J2006||2007[|Hutsi| 
et al.|2009 1 and/or cosmic rays (e.g. |Jasche et al.|2007| l, could affect 



the ionization history at these low redshifts. Overall the uncertainty 
at low redshifts will be dominated by other processes than recom- 
bination physics, and detailed analysis of their relevance is beyond 
the scope of this paper. 



3 PARAMETER ESTIMATION USING COSMOREC 

The difference between ionization history obtained with the orig- 
inal version of Recfast (Seager et al.|l999) and CosmoRec]^] leads 
to differences in the statistics of the CMB that become increasingly 
important at small scales (see Section|2|. Though these differences 
are tiny, in any experiment the smallest cosmic variance limited 
scales are given a huge weight in the likelihood function, and may 
lead to significant biasing when inferring cosmological parameters. 
Our task in this Section is to assess the importance of the deviations 
in recombination calculations on the analysis of future cosmologi- 
cal experiments — our primary approach to this will be to analyse 
simulated data for various experiments using different recombina- 
tion calculations and deduce how well the input model is repro- 
duced. We will compare three different recombination calculations 
CosmoRec, Recfast++ and Recfast++ using a correction function 
(Rubino-Martin et al. 2010), which was obtained using the 'full' 
model given in Table[T]for a cosmology close to the fiducial model 
(given in Table|2]l. We will henceforth call Recfast++ with correc- 
tion function Recfast++cf. 

For this work we rely on CosmoMC (Lewis & Bridle|2002| >, 
a sophisticated Markov-Chain Monte-Carlo (MCMC) package that 
has become the standard code for the analysis of cosmological data. 
At its core is the Boltzmann code CAMB (Lewis et al. 2000 1 which 
calculates cosmological statistics given the proposed set of cosmo- 
logical parameters, and this uses the Recfast v 1.5 system for calcu- 
lating the baryon ionisation fraction and baryon temperature. At the 
simplest level, incorporating CosmoRec into the CosmoMC stack is 
a simple matter of modifying CAMB to replace the Recfast calcu- 
lation. We have made these modifications publicly availabl^] 

In this work we have performed the cosmological analysis 
using several different models, all based around the standard six 
parameter ACDM model. For parameters we use the dimension- 
less baryon and dark matter densities, £l b h 2 and £i c h 2 ; the opti- 
cal depth to last scattering r; the spectral index of the primordial 
power spectrum n s and the logarithm of its amplitude log (10 10 A,) 
(both defined at a pivot scale of k Q = 0.05Mpc~'); finally we will 
use the Hubble parameter Hq (instead of ff). In choosing this pa- 
rameterisation we have enforced a flat universe. We are neglecting 
any secondary contributions to the CMB — including those from 



6 This work was performed using CosmoRec version 1 .3b. 

7 The modified code and instructions can be found at http : //camb . 



the Sunyaev-Zel'dovich effect, gravitational lensing, and inhomo- 
geneous reionisation — as they should be small at the scales we are 
considering (that is / < 2000 in all but one case). 

The additional models we will consider are all extensions to 
ACDM: allowing the primordial Helium fraction Y p to vary; having 
a variable number of neutrino species N v (all of which we assume 
to be massless); and giving the freedom to fit a more complex pri- 
mordial power spectrum by adding the running of the spectral index 
w mn as a parameter. CMB based constraints on Y p and /V,, were re- 
cently obtained by Dunkley et al. (2010) using combined Act and 
Wmap data, demonstrating the power of small scale with full sky 
experiments. In the future such constraints will tighten very much, 
and as we show here the recombination corrections with become 
very important in this case. 

For our simulated data we use an exact realisation of the CMB 
angular powerspectrum calculated for our input model using the 
all_l_exact functionality of CosmoMCQ The noise properties of 
the simulated data are detailed in subsequent sections. We use no 
data (simulated or otherwise) in addition to this CMB data. Our 
input powerspectrum is calculated using the 'full' accuracy setting 
of CosmoRec (see Table [TJ and with all CAMB's accuracy level 
flags set to 4. However when performing the analysis we use the 
'default' CosmoRec accuracy; we give a detailed discussion of the 
implications of this choice in Section]?] 

For analysing cosmic-variance limited data up to a multipole 
of / = 2000 the precision of CosmoMC must be increased to at 
least accuracy_level = 2. This ensures that it is able to repro- 
duce the simulated data at the precision demanded by the cosmic- 
variance limit. If the default accuracy of CosmoMC is used, a log- 
likelihood of 0(1) is obtained for the input model (which should 
have a log-likelihood of zero). In each case discussed we have used 
at least four chains for the MCMC sampling. We have ensured 
tested their convergence using the Gelman-Rubin ( Gelman & Ru- 
|bin|1992] > statistic, in all cases terminating when R — 1 < 0.005. 



3.1 WMAP 

We have performed an analysis of the WMAP 7-year data release 
( |Komatsu et al.|20T0| on its own for the models listed above, using 
the three different recombination calculations. As expected from 
the work of Rubino-Martin et al. (2010), there were no observable 
distinctions between the three, and we confirm that the recombina- 
tion corrections are not important for WMAP. 



3.2 Parameter estimation for Planck 

Planck is expected to be cosmic variance limited up until / ~ 1500 
for temperature observations, (though E-mode polarisation is af- 
fected by instrumental noise at all scales). This increased precision 
means it is much more sensitive to the recombination corrections 
than WMAP ( |Rubino-Martfn et al.|2010l >. 

We have both repeated the analysis of (Rubino-Marti rTet al.| 
|2010[ >, and also extended it by considering models with F p and N v 
included as free parameters. In the simulated data we simply pre- 
sume that the final Planck map is full sky with the same overall 
noise properties as the 143 GHz band (assuming the other frequency 
bands have been used to clean foregrounds). This is the same ap- 
proximation as Rubino-Ma rtin et aT] p010). This means we use a 



inf o/ j r s/cosmorec/ 



See http : //cosmocof fee . info/viewtopic .php?t=231 
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Table 2. The differing recombination calculations could significantly affect the analysis of Planck data. The table below illustrates the biases using each of 
three different recombination codes, in both their absolute deviation and the number of sigmas. Deviations over lcr are highlighted. For the bias in sigmas 
we use multiples of the standard deviation of each distribution, a definition that is robust to their non-gaussianity. Though a more desirable alternative would 
be to turn the confidence limit for each into an effective number of sigmas, this would be highly unreliable as the tails are very sparsely sampled at the large 
deviations observed. 



Parameters Fiducial Recovered CosmoRec Recfast++ Recfast++ w/ correction 

Absolute Sigmas Absolute Sigmas Absolute Sigmas 

ACDM 

Si h h 2 0.0226 0.02260 ± 0.00014 -0.0000 -0.01 -0.0003 -2.10 0.0000 0.01 

Q c h 2 0.112 0.1121 ±0.0012 0.0001 0.05 0.0006 0.48 0.0001 0.04 

H 70 69.98 ±0.60 -0.0200 -0.03 -0.4800 -0.78 -0.0200 -0.03 

r 0.09 0.0905 ± 0.0045 0.0005 0.10 -0.0017 -0.40 0.0003 0.07 

n s 0.96 0.9598 ± 0.0036 -0.0002 -0.06 -0.0120 -3.35 -0.0001 -0.03 

log(10 10 AJ 3.0445 3.0456 ± 0.0089 0.0011 0.12 -0.0091 -1.05 0.0010 0.11 



ACDM + He 

£Vr 0.0226 0.02260 ± 0.00020 0.0000 0.02 2 x 10~ 5 0.10 -1 x 10" 5 -0.04 

a c h 2 0.112 0.1120±0.0013 -0.0000 -0.03 -0.0002 -0.17 0.0001 0.05 

H 70 70.03 ±0.72 0.0300 0.05 0.3800 0.51 -0.0300 -0.04 

r 0.09 0.0905 ± 0.0045 0.0005 0.10 0.0001 0.01 0.0003 0.06 

n s 0.96 0.9601 ± 0.0069 0.0001 0.02 0.0004 0.05 -0.0003 -0.04 

log(10 10 A s ) 3.0445 3.0454 ± 0.0095 0.0009 0.10 -0.0017 -0.18 0.0008 0.08 

Y p 0.24 0.240 ±0.011 0.0000 0.01 0.0230 2.12 -0.0010 -0.05 



ACDM + Neutrinos 

£Vr 0.0226 0.02262 ± 0.00021 2 X 10~ 5 0.08 -0.0000 -0.02 1 X 10~ 5 0.04 

Sl c h 2 0.112 0.1123 ±0.0026 0.0003 0.11 0.0044 1.63 0.0001 0.05 

H 70 70.2 ± 1.5 0.2000 0.12 2.0000 1.21 0.1000 0.05 

r 0.09 0.0907 ± 0.0047 0.0007 0.14 -0.0001 -0.03 0.0005 0.10 

n s 0.96 0.9609 ±0.0081 0.0009 0.11 0.0006 0.07 0.0004 0.05 

log(10 10 A. t ) 3.0445 3.046 ±0.012 0.0020 0.16 0.0040 0.33 0.0010 0.10 

N v 3.046 3.07 ±0.19 0.0200 0.13 0.3300 1.62 0.0100 0.05 



ACDM + Neutrinos + He 

£V; 2 0.0226 0.02261 ± 0.00023 1 X 10~ 5 0.05 5 X 10~ 5 0.24 1 X 10~ 5 0.06 

Q c h 2 0.112 0.1126±0.0038 0.0006 0.17 0.0009 0.24 0.0005 0.13 

H 70 70.3 ± 1.8 0.3000 0.18 0.9000 0.52 0.3000 0.14 

r 0.09 0.0905 ± 0.0046 0.0005 0.11 0.0001 0.03 0.0004 0.10 

n s 0.96 0.9608 ± 0.0083 0.0008 0.10 0.0019 0.24 0.0007 0.08 

log(10 10 A. t ) 3.0445 3.046 ±0.012 0.0020 0.16 0.0000 0.03 0.0020 0.13 

N v 3.046 3.10 ±0.27 0.0500 0.19 0.0900 0.33 0.0400 0.15 

Y p 0.24 0.238 ±0.015 -0.0020 -0.13 0.0190 1.26 -0.0020 -0.11 



ACDM + Running 

Si b h 2 0.0226 0.02260 ± 0.00015 0.0000 0.03 -0.0002 -1.40 0.0000 0.01 

Q c h 2 0.112 0.1 119 ± 0.0013 -0.0001 -0.05 0.0006 0.45 0.0000 0.00 

H 70 70.04 ±0.62 0.0400 0.06 -0.4100 -0.65 0.0100 0.01 

r 0.09 0.0905 ± 0.0046 0.0005 0.11 -0.0000 -0.01 0.0004 0.08 

n s 0.96 0.9601 ± 0.0036 0.0001 0.02 -0.0121 -3.42 0.0000 0.01 

log(10 10 A. t ) 3.0445 3.0454 ± 0.0096 0.0009 0.09 -0.0043 -0.44 0.0009 0.10 

n ran -0.0000 ± 0.0049 -0.0000 -0.01 -0.0058 -1.17 0.0001 0.01 



beam scale of beam = 7.1', and noise of cr 2 ,^,,;,,,, = 1.53x 10 4 /jK 2 
and cr 2 fi bcam = 5.59 x l(T 4 pK 2 . 

The results of our simulations are summarised in Table|2] For 
each parameter in the models listed we give its fiducial value, the 
value constrained by CosmoRec with it standard error, followed by 
the absolute and relative bias for each of the three recombination 
calculations. We now discuss specific cases in more detail. 

3.2. 1 ACDM and ACDM with running 

In the standard ACDM case we reproduce the results of Rubino- 
|Martfn et al.| j2010) finding significant biases in the parameters 



Q. b h 2 and n s when recovered by Recfast++ (~ 2<r and 3<r respec- 
tively). From Fig.[2]we can see that the main effect of the recombi- 
nation corrections on the C/'s is to lower the damping tail further, 
in light of this it is understandable that the bias obtained shifts n s 
in order to lower the small scale multipoles. We note that the biases 
reported here include the total correction to the ionization history 
with respect to the original version of Recfast (Seage r~et al.|1999] >. 
In |Rub ino-Martm et al. (2010), part of the computations were per- 
formed using Recfast v 1.4.2, which allows to account for part of 
the corrections to helium recombination (Wong et al. 2008), using 
fudge factors. Importantly, in comparison to the original version of 
Recfast we find significant shifts for the ACDM + running model 



© 0000 RAS, MNRAS 000, 000-000 



Precise cosmological parameter estimation using CosmoRec 7 



in the parameters Sl),h 2 , n s and n nm (shifts of at least lcr, 3a and la 
respectively), that were not reported in Rubino-M artm et al.| ( |2010[ >. 

3.2.2 ACDM with varying Y p and jV„ 

In the ACDM + Helium model we have allowed variation in the 
primordial Helium fraction Y p rather than simply relying on the 
results from Big Bang Nucleosynthesis (see |Cyburt et al.||2003| 
Steigman 2009 for recent constraints and discussion). Running pa- 
rameter estimations with the Recfast++ reveals a large discrepancy 
in Y v compared to the fiducial model, recovering a value over 2a 
higher. As increasing the Helium fraction removes a larger propor- 
tion of electrons prior to hydrogen recombination, this change can 
partially mimic the lower electron fraction found in a full calcula- 
tion (see Fig. [T|. In fact, virtually all other biases disappear when 
freely varying Y p . With ~ 0.5<x only H is significantly biased in 
addition for this case. 

The CMB is sensitive to the effective number of neutrino 
species N v primarily because they affect the background radiation 
density of the Universe. Any significant deviation from the theoret- 
ical value of 3.046 (Mangano et al. 2005) could indicate a further 
generation of neutrinos, or potentially some other source of back- 
ground radiation. As illustrated in Table [2] to infer this value cor- 
rectly requires an accurate calculation of the recombination history. 
To do otherwise can lead to a bias of up to 1.6cr in the recovered 
value (as shown by using Recfast++), as well as a similar shift in 
Sl),h 2 , and a slightly smaller change in H . The high dimensionality 
of the parameter space, and the numerous shifts makes it difficult 
to attribute them to any definite physical effect. We also allow both 
Y p and N v as free parameters, in an eight parameter model. This 
results in a smaller set of shifts, with a higher value of Y p acquiring 
the most significant shift (at 1.3cr). 



nation corrections for any particular experiment, or combination of 
experiments, we will use a simulated, full-sky, experiment cosmic 
variance limited until some maximum multipole used in our analy- 
sis, usually Z = 2000. Apart from this we perform a similar analysis 
to the one of the previous Section. Our results are summarised in 
Table |3] 

Overall the obtained set of biases for each combination of pa- 
rameters are similar to the Planck cases shown previously. How- 
ever, the relative shifts are much greater in every case. For instance 
in the plain ACDM model, the bias increases to over 5 and la for 
Qi,h 2 and n s respectively (also see Fig. [5}. In this case the signif- 
icant driver of this increase is the overall tightening of the con- 
straints from the added data. However, in many cases there is also 
an absolute shifting of the biased parameters. This happens as the 
weight afforded to the largest multipoles eliminates the previously 
favourable regions, shifting them further away. This is most obvi- 
ous in the shift of Ho in the ACDM + Helium model (see Fig.|5]as 
well as Table [3j. 

In Tabled we also include the results from the ACDM model 
using simulated data from an experiment cosmic variance limited 
up to / = 3000. In this case the uncertainty on Qkh 2 reduces by 
around a half, which increases the associated bias to over lOcr. For 
n s , however, the distribution width remains comparable, but an in- 
crease in the absolute shift, moves the bias to over 12cr. 

Here it is important to mention that at / ~ 3000 the additional 
ambiguities caused by SZ clusters, are expected to mainly impact 
the TT power spectra, while leaving the EE power spectra less con- 
taminated. This holds the potential that foreground subtraction and 
the removal of secondary anisotropies up to / ~ 3000 could become 
feasible in the future, in particular for the £-mode power spectrum. 
As our results clearly show, in such cases the recombination cor- 
rections will be extremely important. 



3.2.3 Parameter estimation using Recfast++cf 

As the results in Table [2] indicate, at the precision of this Planck 
simulation there is no significant difference between the recovered 
result for either of the CosmoRec or Recfast++cf cases for any cos- 
mological model. This confirms that a correction function approach 
in principle is sufficient for the analysis of Planck data, as already 
pointed out by Rubino-Ma rtm et al.| l|2010). However, given that 
one run of CosmoRec with 'default' setting takes ~ 1.3 seconds or 
about ~ 15% of the total CPU time of CAMB per modeQ it has 
now become possible to explicitly run the full recombination code, 
without large penalty. 

3.3 Cosmic variance limited parameter estimation 

Future experiments such as SPTpol and ACTPol will be much 
more sensitive than Planck at small scales. This means that the 
combined Planck plus SPTpol and ACTPol data sets will be sam- 
ple variance dominated to much smaller scales than currently ac- 
cessible^] Rather than exploring the consequences of the recombi- 

9 This is with a timing of 7.5 seconds per CAMB call (at accuracy.level 
= 2) on a single processor. Using OpenMP, CAMB can be parallelised to 
reduce its overall time, however, CosmoRec does not benefit from paralleli- 
sation at the moment. 

10 Up to / ~ 2000 for £-mode polarisation in the case of Planck plus the 
ACTPol wide survey. ACTPol deep may go as far as / ~ 4000. See |Niemack| 
[eTaT|(20T0| 



4 EFFECT OF PRECISION IN THE RECOMBINATION 
PHYSICS 

In Section [3] we have used CosmoRec with the 'default' setting to 
perform model-by-model computations of the recombination pro- 
cess. However, it is clear that this leads to some residual errors 
with respect to the most precise recombination calculation (model 
'extreme' in Table[T]). In terms of the CMB power spectra this 'ap- 
proximation' leads to very small uncertainties, as discussed already 
in Section|2](see Fig.|2]for details). Here we ask at what level these 
differences actually matter for future CMB parameter estimation, 
using explicit runs of CosmoRec. 



4.1 Changes to the freeze-out tail of recombination 

In Section l2T2l we discussed the effect of the number of shells calcu- 
lated on the low redshift ionisation fraction. We expect the relative 
insensitivity of the CMB power spectrum, means that the number 
of shells included is not of great importance. Using CosmoRec this 
can be explicitly verified by varying the number of hydrogen lev- 
els included into the computation of the effective rate coefficients. 
In Fig. [6] we can see that even in the cosmic variance limited case 
there is no significant bias introduced when reducing the number 
of shells down to 100. For Planck we expect 100 shells to be suffi- 
cient, confirming the analysis of Rubino-M artm et al.| ( |2010[ >. 
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Table 3. Deviations between the different recombination calculations for an experiment cosmic-variance limited until / = 2000 (except in the final case listed). 
The relative deviations are significantly larger than those in a Planck like experiment (Table|5J. Any deviations over 2tr have been highlighted. 



Parameters 


Fiducial 


Recovered 


CosmoRec 


RECFAST+ + 


Recfast++ w/ correction 








Absolute 


Sigmas 


Absolute 


Sigmas 


Absolute 


Sigmas 


ACDM 


















n b h 2 


0.0226 


0.022594 ± 0.000039 


-6 x ur 6 


-0.15 


-0.0002 


-5.35 


-1 X 10~ 6 


-0.01 


n c h 2 


0.112 


0.11201 ±0.00051 


1 X 10~ 5 


0.01 


-0.0004 


-0.77 


-2 X 10~ 5 


-0.04 


Ho 


70 


69.99 ± 0.22 


-0.0100 


-0.03 


0.0000 


0.00 


0.0100 


0.05 


T 


0.09 


0.0901 ± 0.0020 


0.0001 


0.03 


-0.0013 


-0.64 


0.0000 


0.02 


n s 


0.96 


0.9598 ± 0.0016 


-0.0002 


-0.14 


0.0122 


-7.48 


0.0000 


0.01 


log(10 10 A s ) 


3.0445 


3.0446 ± 0.0042 


0.0001 


0.04 


-0.0071 


-1.72 


0.0000 


0.01 


ACDM + He 


















n b h 2 


0.0226 


0.022602 ± 0.000059 


2 X 10~ 6 


0.03 


0.0001 


1.78 


1 X 10~ 6 


0.02 


a c h 2 


0.112 


0.11201 ±0.00051 


1 x 10~ 5 


0.02 


-0.0006 


-1.06 


6 x 10~ 5 


0.11 


Ho 


70 


70.01 ± 0.23 


0.0100 


0.03 


0.6000 


2.50 


-0.0100 


-0.06 


T 


0.09 


0.0901 ± 0.0021 


0.0001 


0.04 


0.0002 


0.08 


0.0001 


0.04 


»s 


0.96 


0.9601 ± 0.0025 


0.0001 


0.03 


0.0023 


0.87 


-0.0000 


-0.00 


log(10 10 A s ) 


3.0445 


3.0449 ± 0.0043 


0.0004 


0.10 


-0.0005 


-0.10 


0.0006 


0.13 




0.24 


0.2405 ± 0.0035 


0.0005 


0.13 


0.0250 


7.24 


0.0002 


0.07 



ACDM + Neutrinos 



il b h 2 


0.0226 


0.022597 ± 0.000058 


-3 x ur 6 


-0.05 


5 x ur 5 


0.93 


-3x 10- 6 


-0.06 


n c h 2 


0.112 


0.11203 ±0.00099 


3 x 10~ 5 


0.03 


0.0045 


4.77 


-0.0000 


-0.00 


Ho 


70 


69.99 ± 0.45 


-0.0100 


-0.02 


2.2800 


5.40 


-0.0600 


-0.14 


T 


0.09 


0.0902 ± 0.0021 


0.0002 


0.08 


-0.0002 


-0.11 


0.0000 


0.02 


n s 


0.96 


0.9600 ± 0.0027 


0.0000 


0.01 


0.0007 


0.26 


-0.0002 


-0.08 


logao 10 ^) 


3.0445 


3.0449 ± 0.0049 


0.0004 


0.08 


0.0066 


1.37 


0.0002 


0.03 


Ny 


3.046 


3.046 ± 0.058 


0.0000 


0.01 


0.3460 


6.35 


-0.0050 


-0.08 


ACDM + Neutrinos + He 
















Q b h 2 


0.0226 


0.022602 ± 0.000061 


2 X 10- 6 


0.04 


0.0001 


1.94 


2 x 10~ 6 


0.03 


a c h 2 


0.112 


0.1121 ±0.0014 


0.0001 


0.09 


0.0007 


0.50 


0.0001 


0.06 


Ho 


70 


70.04 ± 0.54 


0.0400 


0.07 


1.1000 


1.94 


0.0100 


0.01 


T 


0.09 


0.0900 ± 0.0022 


0.0000 


0.02 


0.0001 


0.04 


0.0001 


0.04 


n s 


0.96 


0.9601 ± 0.0028 


0.0001 


0.03 


0.0032 


1.11 


-0.0000 


-0.01 


log(10 10 A.,) 


3.0445 


3.0448 ±0.0051 


0.0003 


0.05 


0.0016 


0.31 


0.0003 


0.07 


Ny 


3.046 


3.054 ± 0.087 


0.0080 


0.09 


0.0890 


1.01 


0.0040 


0.04 


Y P 


0.24 


0.2396 ± 0.0053 


-0.0004 


-0.07 


0.0206 


3.99 


-0.0005 


-0.09 



ACDM + Running 



a b h 2 


0.0226 


0.022594 ± 0.000043 


-6 x 10~ 6 


-0.15 


0.0002 


4.05 


-2 x 10~ 6 


-0.05 


Cl c h 2 


0.112 


0.11203 ± 0.00055 


3 x 10~ 5 


0.06 


-7 x 10~ 5 


-0.13 


4x 10~ 5 


0.07 


Ho 


70 


69.99 ± 0.23 


-0.0100 


-0.06 


-0.1000 


-0.42 


-0.0100 


-0.05 


T 


0.09 


0.0901 ± 0.0022 


0.0001 


0.04 


0.0006 


0.28 


0.0001 


0.03 


n s 


0.96 


0.9598 ± 0.0017 


-0.0002 


-0.10 


0.0118 


6.98 


-0.0000 


-0.00 


log(10 10 A. t ) 


3.0445 


3.0449 ± 0.0047 


0.0004 


0.08 


-0.0000 


-0.00 


0.0004 


0.09 


'^run 





-0.0001 ± 0.0027 


-0.0001 


-0.05 


-0.0088 


-3.41 


-0.0001 


-0.03 



ACDM (CV limited up to / max = 3000) 



a b h 2 


0.0226 


0.022600 ± 0.000019 


-0.0000 


-0.03 


-0.0002 


10.32 


3 x 10~ 6 


0.18 


Cl c h 2 


0.112 


0.11202 ±0.00047 


2 x 10~ 5 


0.03 


-0.0002 


-0.38 


4x 10~ 5 


0.09 


Ho 


70 


69.99 ±0.19 


-0.0100 


-0.03 


-0.1000 


-0.52 


-0.0100 


-0.06 


T 


0.09 


0.0901 ± 0.0020 


0.0001 


0.07 


-0.0021 


-1.18 


0.0000 


0.01 


n s 


0.96 


0.9600 ± 0.0014 


0.0000 


0.03 


-0.0175 


-12.66 


0.0001 


0.08 


log(10 10 A s ) 


3.0445 


3.0448 ± 0.0042 


0.0003 


0.08 


-0.0060 


-1.62 


0.0003 


0.06 



4.2 CosmoRec 'default' versus 'full' setting 

As remarked upon in the previous section, our simulated data has 
been generated with the 'full' setting of CosmoRec, which takes 
around 50 seconds to evaluate the recombination history. However, 
this is too slow for use during an MCMC analysis, as it would 
make one execution of CAMB even with highest accuracy setting 



about seven times slower. By using the 'default' setting we are in- 
troducing differences in the power spectrum of around 0.05% (at 
/ ~ 2000; cf. Fig. E), though small, this induces slight changes to 
likelihood surface, shifting the preferred parameters. This means 
that even the CosmoRec case exhibits slight biases when looking 
at cosmic variance limited data (see Table [3}. However, these are 
limited to around 0. 15cr. This deviation can be reduced by running 
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Figure 4. Recovered constraints from a simulated experiment, cosmic variance limited in both temperature and polarisation until Z max = 2000. This plot shows 
the constraints obtained for the six parameter model, with the dashed lines and crosses marking the input values in the one- and two-dimensional plots. The 
contours mark the 68% and 95% confidence limits of each distribution. The three cases illustrated show that using both CosmoRec and corrected Recfast++ 
produce nearly identical probability distributions which faithfully recover the input model. In contrast, the standard Recfast++ model significantly biases the 
parameters (most notably Cl\,h 2 and n s ). 



CosmoRec at higher accuracy. Using the 'HyRec' setting (which 
does not include the Hei diffusion correction, but takes 4s/4d-ls 
two-photon decay and 3s/3d-ls Raman events into account; see Ta- 
ble[T]for more details) the maximum bias is reduced down to 0.09<r, 
and if the He i diffusion corrections is also included this decreases 
to 0.06tr. In both cases the runtime of CosmoRec is increased only 
by about 1 second, a relatively small penalty. We conclude that for 
Planck parameter analysis the CosmoRec 'default' setting is suf- 
ficient, however, if needed additional corrections can be included 
without compromising the computation time very much. Further- 
more, our estimation shows that the He i diffusion correction leads 
to a small additional modification, which can be safely neglected. 



4.3 Iterative use of Recfast++cf 

One alternative to using an explicit CosmoRec calculation (or any 
other more detailed recombination code) is to resort to an incre- 
mental approach, based on Recfast++ with correction function. In 
the absence of any prior knowledge about the appropriate cosmol- 
ogy, this would commence with an initial analysis of the data using 



Recfasth — h . Then the best fit parameters are used to generate a cor- 
rection function for that cosmology using the most precise recom- 
bination calculation (e.g. Chluba & Thomas 2010, Ali-Hai'moud & 
|Hirata|2010a[ >, which should be correct within a nearby region in 
parameter space. The data is then re-analysed with Recfast++ us- 
ing this correction function. This process can be iterated until the 
results are convergent. 

We have shown the results of this process in Fig. [7] where we 
illustrate the Q. b h 2 -n s plane (projected from the full ACDM param- 
eter space), using the cosmic -variance, / raax = 2000 simulated data. 
The initial results from Recfast++ are highly biased, as discussed 
in the previous section. However the first iteration significantly cor- 
rects this, with the largest residual being the 0.26<x shift in Cl/Ji 2 , 
indicating that 'corrections-to-corrections' are very small already. 
A further iteration reduces this shift down to only 0.03cr, suggest- 
ing that in this cases two iterations are sufficient. This could likely 
be reduced down to only one iteration if our cosmological knowl- 
edge is sufficient to calculate a suitable initial correction function, 
eliminating the Recfast++ only step. 

Furthermore, at the initial step one could use Recfast vl.5 that 
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n h h 2 
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Figure 5. Constraints with our cosmic-variance limited experiment for / < 2000, for the seven-parameter model of ACDM plus an unknown primordial Helium 
content, F p . The extra degree of freedom shifts the bias of the standard Recfast++ calculation into a dramatically higher value of Y p . 



O b h 2 




— 500 Shells 

— 300 Shells 

— 100 Shells 



0.02250 0.02265 0.022S0 

n s 




0.02250 0.02265 0.02280 0.954 0.960 0.966 



Figure 6. The projection of the Hj } h 2 -n s plane from the 6-parameter ACDM 
model, analysed using a maximum of 100, 300 and 500 shells in CosmoRec. 
As can be seen for an experiment cosmic-variance limited at / = 2000 or 
less, using a maximum of 100 shells is sufficient. 



mimics some of the recombination corrections using fudge func- 
tions that were calibrated on the results of earlier precise recombi- 
nation calculations (Rubino-Martm et al. 20101. However, since the 
fudge function approach ( Wong et al. 2008) in essence is equivalent 
to a correction function approach ( [Rubino-Martm et al.|2010| t, one 
does not expect any additional improvement in this case. Once we 
are unwilling to run improved recombination codes (i.e. CosmoRec 
or HyRec), a correction function approach will lead to a very simi- 
lar approximation at re-calibrated fudge functions to Recfast. 



5 CONCLUSIONS 

In this work we have demonstrated in detail how previously ne- 
glected physical processes during the recombination epoch affect 
the analysis of precise CMB data from Planck and beyond. We 
have successfully incorporated CosmoRec into the widely used pa- 
rameter estimation code CosmoMC and the Boltzmann code Camb. 
With this we were able to show the parameter biases introduced in 
a comprehensive set of models when neglecting these corrections. 

We confirm the significant biases in the standard six-parameter 
ACDM model previously reported by Rubino-M artm et aL"|(2010) 
for case of Planck . In addition to considering a wider range of 
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0.02250 0.02265 





0.02235 0.02250 0.02265 



Figure 7. Successive iterations of the Recfast++ correction function in 
the Qi,/i 2 -«j plane (taken from a full six-parameter ACDM analysis). The 
'No Correction' plot, calculated using Recfast++, exhibits significant bi- 
ases, however, from this we can calculate the first correction function. This 
correction function can be used in Recfast++, to produce improved con- 
straints. This iteration is much less biased, with the largest change exhib- 
ited a 0.26cr shift in fl/,/1 2 . A further iteration with a recalculated correction 
function, removes the last remaining biases. 



cosmological models, we also extend this investigation to a future 
dataset (comparable to Planck plus ACTPol/SPTpol) using full 
model-by-model computations of the ionization history with Cos- 
moRec (see Tables [2] and |3] for details). There are several biases of 
particular note, which have not been mentioned elsewhere: 

• Allowing a running of the spectral index yields a bias of over 
lcr in n mn for Planck. For an experiment that is noise-free up to 
/ = 2000 this increases to over 3rx. 

• If the primordial Helium fraction is to be inferred using CMB 
data alone it may be biased by up lcr. If uncorrected this would er- 
roneously indicate a significant tension with predictions from Big- 
Bang Nucleosynthesis. 

• Adding N v as a free parameter would bias its estimated value 
high by up to \.6cr for Planck , however, when also varying the 
helium fraction this decreases to 0.3cr only. 

Our analysis indicates that for the current precision of com- 
bined WMAP plus ACT data, the bias in Y p corresponds to about 
0.5cr. This is small compared to the ~ 2<r tension with the BBN 
value reported by Dunkley et al. (2010). However, in the near fu- 
ture the error in the measurement should decrease and the corre- 
sponding bias in Y p caused by a neglect of detailed recombination 
physics should exceed the level of lcr, even prior to the final release 
of Planck data. 

Although small scale CMB data is able to internally break 
many of the degeneracies of current cosmological analyses, it will 
still be desirable to supplement future CMB data with complemen- 
tary probes such as Type la supernovae data {A manullah et al.| 
|2010| > or BAO measurements (Schle gel et al.|2009) . We have not 
included any additional data in this work, and whilst it may serve 
to reduce the overall biases obtained (if using a Recfast analysis), 
the tension between datasets is likely to artificially tighten parame- 
ter constraints. 

In order to focus on the importance of the recombination 



corrections we have necessarily neglected many instrumental and 
physical systematics: 

• Both the SZ effect and gravitational lensing will significantly 
add to the small scale CMB. We have neglected these in our analy- 
sis assuming that for / < 2000 they will be negligible or sufficiently 
well understood. 

• Our work uses a simple reionisation model which is as- 
sumed to be known (other than a single parameter giving the opti- 
cal depth). A more complicated homogeneous reionisation history 
should only affect the large scale multipoles, whereas any realis- 
tic inhomogeneity will not significantly effect multipoles / ~ 2000 
(see e.g., |Zahn et al.|2005^ 

• We have assumed a simple, perfectly known, noise model in 
this work. In reality beam, gain and other uncertainties will have to 
be considered (see e.g. Col ombo et al.|2009| >. 

The 'default' accuracy setting for CosmoRec is designed to be 
a fast enough to be used for MCMC analysis of cosmological data, 
and this work shows that it is sufficiently accurate for the unbiased 
analysis of Planck data. We also discuss an iterative approach us- 
ing Recfast++ with a correction function. Though this approach is 
accurate (see Sect. |4.3fr , it is simpler and quicker overall to directly 
use the full CosmoRec analysis. For the analysis of future Planck 
plus ACTPol/SPTpol data, we show that CosmoRec with 'default' 
setting is unbiased at the 0.1 5<x level, which can be significantly 
reduced by running at a higher accuracy setting (leading to a small 
runtime increase). 

It is important to note that a final cross-validation of the Cos- 



moRec outputs with independent recombination codes I Switzer & 



|Hirata|20081 |Hirata|2008||Grin & Hirata 20T0l |Ali-Haimoud et al. 
|2010[ > will be very important, for confirming their veracity. Never- 
theless, initial comparisons mean that we do not expect any signif- 
icant changes (at least for the experiments simulated in this work). 

We would also like to mention that in this work we have 
demonstrated the importance of the total corrections to the ioniza- 
tion history with respect to the original Recfast code jSeager et al.| 
1999). This correction was obtained in a common effort by several 
independent groups, and part of these corrections are now included 
by the latest version of Recfast using calibrated correction func- 
tions, based on the more detailed recombination computations. As 
we demonstrated here, such an approach is similar to an iterative 
scheme using a precomputed correction function for the original 
Recfast computation (see Sect. |4.3[ >. However, we also argued that 
for runtimes achieved with recent recombination codes, such an ap- 
proximate approach can be avoided. 
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